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Using the method of uniqueness a two-loop massless propagator Feynman diagram with a nonin- 
teger index on the central line is evaluated in a very transparent way. The result is applied to the 
computation of the two-loop polarization operator in reduced quantum electrodynamics. 



Introduction - The exact analytical computation of 
multi-loop Feynman diagrams is of crucial importance for 
the evaluation of renormalization group functions, i.e., j3 
functions and anomalous dimensions of fields. Since the 
early days of quantum field theory a variety of meth- 
ods have been developed, and often combined, in order 
to achie ve t his task, e.g., the G egen bauer polynomial 
techniqu^l^, integration by part^^lil^ ^^^^ ^j^g niethod 
of uniquenesP'SHZI The latter allows, in principle, the 
computation of complicated Feynman diagrams using se- 
quences of simple transformations. A given diagram 
is straightforwardly integrated once the appropriate se- 
quence is found. The task of finding such a sequence for 
a given diagram is, however, nontrivial, see Ref. |H| for a 
review. 

One of the basic building blocks of multiloop calcula- 
tions is the two-loop massless propagator diagram: 

J(ai,Q!2,a3,a4,a5) = (1) 
f f d^fci d^fca 

with arbitrary indices and external momentum p in 
a Euclidean space-time of dimensionality D, Fig. [T] see 
Ref. |9j for a historical review on this diagram. When 
all indices are integers the diagram of Eq. is well 
known and easily calculated. Its evaluation for arbi- 
trary indices is however highly nontrivial: the results can 
be represented^ as a combination of twofold series. In 
some particu lar cases, however, the results can be ob- 

tainedaMnmni 

in significantly simpler form. In Ref. 
a class of complicated diagrams, with two integer indices 
on adjacent lines and three other arbitrary indices, has 
been computed exactly on the basis of a new development 
of the GP technique. The latter approach is technically 
involved and the result is expressed in terms of a general- 
ized hypergeometric function, 3^2 with argument 1. For 
this class of diagrams, similar results have been found in 
Ref. [13] using an ansatz to solve the recurrence relations 
for the two-loop diagram arising from IBP. In this brief 
report we consider the simplest but important diagram 
belonging to this class, whose coefficient function reads 

„2(2-A) jj 

IiX)^^^^J{hl,l,l,X), A=--l, (2) 




FIG. 1: Two-loop massless propagator diagram. 
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FIG. 2; Two-loop vacuum polarization diagrams. 

where the index 015 has been restricted to A, all other in- 
dices being 1. This diagram has been already calculatecP 
(see also discussions in Ref. [15^) but it seems that our 
evaluation is simpler and more transparent. As will be 
seen in the following, the diagram can indeed be straight- 
forwardly integrated using the method of uniqueness in 
momentum space with the help of a simple but inge- 
nious three-step transformation. As an application, we 
compute the two-loop polarization operator (Fig. [2| in 
reduced quantum electrodynamicl^S. 

Definitions and notations - In what follows we use di- 
mensional regularization. All calculations are performed 
in a Euclidean space-time of dimensionality D = 2 + 2X 
which may either be even dimensional (A — 1) or odd 
dimensional (A — > 1/2). Diagrams will be analyzed in 
momentum space. For a given diagram, the integrations 
are over loop momenta and the lines are simple power 
laws of the form: 1/A;^" where a is the index of the line. 
The index of a diagram is defined as the sum of the in- 
dices of its constituent lines. A line with an arbitrary 
index a can be represented graphically as: 




In order to simplify notations, arrows will not be dis- 
played in the following but they are implied, e.g., Fig. [T] 
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With these definitions and graphical notations, chains reduce to the product of propagators: 

a + P 



a P 
• 



(4) 



On the other hand simple loops involve an integration: 



/3 




a + 13- D/2 



A! m a{c.)a{p) T{Dl2-a) 

a[a + p — 0/2) i (a) 



In momentum space a vertex and a triangle are said to be unique if their indices are equal to Z?/2 = 1 + A and 
D = 2 + 2A, respectively. They are related to each other by the uniqueness (or star-triangle) relation: 




^2 



(6) 



where — D /2 — ai is the index dual to ai. Finally, for an arbitrary triangle (unique or not) the following recurrence 
relation is obtained from integration by part 4^'^ 



{D - a2- as. - 2a:-,) 




a2 



a 3 



(7) 



where ± on the right-hand side of the equation denotes the increase or decrease of a line index by 1 with respect to 
its value on the left-hand side. In the following, in order to simply notations, we will assume that lines with no index 
are ordinary lines. In momentum space ordinary lines have index a = \. 

Calculation of the diagram - With the help of the above notations and identities we proceed on calculating /(A). 
The first transformation consists in replacing the central line by a loop^'^, Eq. ([5|, in order to make the right triangle 
unique. The uniqueness relation, Eq. ([6]), can then be used. In graphical notations this reads 



"(^'^'^'1'^) - - .W^i(l,2A) 




p 



,2(1-A) 



(8) 



Finally, using integration by parts, Eq. ([t]), the last diagram is reduced to sequences of chains and simple loops which 
can immediately be integrated with the help of Eqs. Q and Q: 

A + (5 

^>^ + 5 r ^ — ^ — ^.-t-^A + (5 



(-2(5) h/ V = 2(A + 5) 

TT°2{\ + 5 ) 
~ p2(l+2<5) 




A + 5+1 



A + (5 + l 



^(1,1) A{\ + 5 + l,\ + 5)- A{\ + 5 + l,l + 6) 



where the parameter 5 has been introduced. The function ^(1, 1) and the bracketed terms in Eq. (9b) read 

r2(A)r(l- A) 



A(l,l) 
(A + <5) 



r(2A) 

r(-^) r(A - <5)r(i 2(5) 



T{\ + 5) T{l + 5)T{\-25) 



r(l - (5)r(l + (5)r(A + 25)T{X - 25) _ 
r(l - 2(5)r(l + 26)T{\ + (5)r(A - 5) ^ 



(9a) 

(9b) 

(10a) 
(10b) 



J 2 
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At this point it is convenient to use the following product expansion of the Gamma function: 



r(a; + e) = r(a;) exp [ ^^^'^''^Ha 



fc=i 



where i/''''^'' is the polygamma function of order k. From Eq. the following relation is obtained 



oo 

T{x + e)T{x-e)^r^{x) exp [ 2 ^ V'^^^^^a^) 



-2m 



(2to)! 



Making use of Eq. ( 12 ) in the bracket of Eq. ( 10b ) yields 



oo 

exp [ 2 ^ (22" - 1) [v;(2'"-i)(A) - 



r2m 



(2to)! 



3<5^ 



V''(A)-^'(1) +0{6^) 



(11) 



(12) 



(13) 



Substituting back Eq. ( 13 ) in ( 10 ) and performing the remaining S expansion yields 

A 

_ r(A)r(i-A) 



(A + <5) 





-3d 








• 


1 - r(A) 


• 


3 EqT^|9| \^ 


) 



r(2A) 



V''(A)-^'(l) 



(14) 



A 



where, in the last step, Eq. ([9| has been used and S sent to zero. Substituting the final result of Eq. (14) in Eq. ([s]) 
and using Eq. ([2|, we obtain the advertised resullPE^for the coefficient function: 



/(A) =3 



r(A)r(i - A) 
r(2A) 



^'(A)-v/(i) 



(15) 



where ip'ix) is the trigamma function. In the even-dimensional case (A — >■ 1 or Z? — >■ 4) the well-known result 
/(I) = 6C(3), is obtained. On the other hand, in the odd-dimensional case (A 1/2 or D — >■ 3), which is one of the 
cases of interest to Ref. |12|16j . the result reads /(1/2) = 6tt({2). 

Application - We now focus on the computation of radiative corrections to the polarization operator 11'^'^ (g) = 
n(g2) {g^'^q^ — q'^q'^) in reduced quantum electrodynamics^^ (RQED^ ^d^). The latter describes the interaction of a 
photon field living in dimensions with a fermion field living in a reduced space-time of c?e dimensions (c?e < d^). 
Within dimensional regularization, the computation of Feynman integrals in such a reduced theory can be carried out 
by introducing two epsilon parameters, e-y and Eg, such that c?^ = 4 — 2e^ and de = 4 — 2ee — 2e^, respectively. In 
Ref. |16| the corrections up to two loops (see Fig. [2] where the corresponding diagrams were displayed) were computed 
for an arbitrary KQEDd^ d^ using the general result of Ref. [2] for /(A). The resulting expression is rather cumbersome. 
Here, we focus on the case of RQED4 . In the limit e-^ — ^ and using Eq. ( 15 ) we obtain the following simpler and 
more explicit formulas (see definitions in Ref. [16]): 

, e4r(A)r2(l + e^) 16r(l + A)r(l-A) 



2IL2aiqn 



(47r 



X <^ A2 — + V 



r(3 -I- 2A) 
3A „ 2 



n2fc(g') = d 



e4r(A)r2(i + £ 



1 + 2A/ 2 1 + A 

16r(l + A)r(l - A) 



(4^)3+A-2s-, (^2)l--^+2'^. 
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X -A^ 



— +ip- 



1 + 2A 



r(3 + 2A) 

A _ 1 
2^2 



3 

2A 



+ O(e^) } , (16a) 



+ ^A(l + A) [^'(A) - ^'(1)] + 0(£^) ,(16b) 



where A = 1 — Ee and ^ — 3ip{2X) — 2tp{X) + 2^{1 — A) — 3tp{l). The one-loop and total two-loop contributions 
therefore read 



Uliq^) = -d 
U2{q^) = d 



e^r{x) r(i + A)r(i-A) 

(4^)i+A (^2)i-A (l + 2A)r(2A) 

e^r(x) i6r(i + A)r(i-A) 



(47r)3+A (q2) 



1-A 



r(3-l-2A) 



CliX); 



^ + ^A(1 + A)[^'(A)-^'(1) 



(17a) 
(17b) 
(17c) 
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From Eqs. (17) we see that ni((7^) and n2(g^) are finite as long as A 7^ 1. We can then replace by 4TTa in (17) 
which yields 



n2(g') = - 



•A(l + A) 



Ci(A)ni(q2), 



(18) 



It should be noted that these results cannot be used for QED4 which can be reached from a general KQEDd^^d^ (see 
Ref. [H]) by first fixing Ee — and then taking the limit e-^ 0. The results ( 17) are singular in the limit A — 1 but 
this limit corresponds to = and Eg 0, which does not lead to QED4. 

The total, up to two loops, gauge-field self-energy in RQED4_d^ (e-y — and arbitrary Sg) may then be written as 



C(A) 



7rA(l + A) 



Ci(A) 



27r 



^'(A + 2)-V^'(l) 
I 



1 + A (1 + A)^ 



(19) 



For A = 1/2, i.e., in the case of RQED4_3 {e^ = and 
Ee = 1/2) which corresponds to an ultrarelativistic model 
of graphen^iSl (a 2-brane), we reproduce the basic result 
of Ref. dg P 



C7i(l/2) 



97r2 - 92 
24 ' 



C{l/2) 



92 - 9tt^ 
IStt 



(20) 



This coefficient is small, C(l/2) sa 0.056, in qualitative 
agreement with some results obtained in the nonrelativis- 
tic limilPSl (see however Refs. [21]) as well as experimental 
result^^ where C(l/2) corresponds to an interaction cor- 
rection coefficient to the optical conductivity of undoped 



graphene. 

As a next step of our future investigations we would 
like to evaluate the fermion self-energy of general 
RQEDjj d^ in the ultrarelativistic limit. 
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In coordinate space, it corresponds to the insertion of a 

point into this line (see the table of such transformations 

in Ref. [3| and also Ref. [S] for a review). 

A coefficient similar to Ci (1/2) has been obtained earlier in 

|17| in the framework of the 1/A'^ expansion of QED3. The 

similarity of the coefficients comes from the corresponding 

similarity of the effective photon propagators in RQED4,3 

(see Eq. (9) in [TB]) and in the 1/A'' expansion of QED3 

(see Eq. (2.16) in [^) at large N. 



